Reduced effect of impurities on the universal pairing scale in the cuprates 
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We consider the effect of non-magnetic impurities on the onset temperature T* for the d— wave 
pairing in spin-fluctuation scenario for the cuprates. We analyze intermediate coupling regime when 
the magnetic correlation length ^/a > 1 and the dimensionless coupling u is 0(1). In the clean limit, 
T* K, 0.02vf /a in this parameter range, and weakly depends on ^ and u. We found numerically that 
this universal pairing scale is also quite robust with respect to impurities: the scattering rate Fcr 
needed to bring T* down to zero is about 4 times larger than in weak coupling, in good quantitative 
agreement with experiments. We provide analytical reasoning for this result. 



Introduction. The issue of the pairing glue in the 
cuprates is still one of the hottest topics in the physics of 
strongly correlated electrons. Many researchers believe 
that the pairing glue is a spin fluctuation exchange, at 
least in the overdoped and optimally doped cuprates. At 
such dopings correlations are relatively strong, but not 
enough so to break a large normal state Fermi surface 
(FS) into the hole and electron pockets. 

Recently, Abanov et al. analyzed the pairing prob- 
lem in the cuprates within spin-fluctuation scenario, un- 
der the assumption that the FS in the normal state is 
largcii Within this scenario, the onset of the pairing at 
T = T* marks the development of the pseudogap phase 
with strong thermal fluctuations of the pairing gap, and 
the true superconductivity emerges at a smaller Tc < T* , 
when thermal fluctuations become weaker. These au- 
thors found a smooth crossover between the limit when 
the interaction U is smaller than Vf/a, and the pair- 
ing is confined to a vicinity of hot spots (points on the 
FS separated by (tt, tt)), and the limit of strong interac- 
tion, when the entire FS is "hot" {vf is the bare Fermi 
velocity as obtained in band theory, and a is the in- 
teratomic spacing, Vf/a ~ leV). The onset temper- 
ature for the pairing, T*, scales as T* ~ {vf/a)u for 
u ~ 3Ua/8TTVf ^ 1, and as T* ^ {vf/a)(l/u) at large 
M, and weakly depends on ^ for ^ > a. For interme- 
diate values of u = 0(1), which are mostly relevant to 
the cuprates, T* goes through a shallow maximum and 
is approximately 0.02vf/a. The same pairing scale was 
obtained in FLEX calculations for the Hubbard model^, 
in two-particle self-consistent calculations^, in dynami- 
cal cluster approximation^, and in cluster DMFT— . The 
good agreement between all these results is strong in- 
dication that T* ~ 0.02w//a ~ 200 - 250A" is indeed 
the universal pairing scale in optimally doped cuprates. 
At smaller dopings, this scenario breaks down because of 
electron localization which gives rise to precursors to hole 
and electron FS pockets already in the normal state. The 
pseudogap temperature T* then becomes a scale at which 
the system develops such precursors, while the pairing 
emerges at a smaller temperature due to interaction be- 
tween electron pockets. 

The subject of this communication is the analysis 



of how the universal pairing scale T* in optimally- 
doped cuprates is affected by non-magnetic impurities, 
which are pair-breaking for unconventional superconduc- 
tors. In near-optimally doped cuprates, concentrations 
of dopants are quite substantial, and potential random 
scattering off dopants could significantly reduce T*. At 
weak-coupling, which in our case corresponds to small 
u and ^ a, non-magnetic impurities in a c?— wave 
superconductor suppress T* in the same way as mag- 
netic impurities in a BCS superconductor, and T* is 
given by Abrikosov-Gorkov (AG) formula- logTQ*/r* = 
^-(1/2 + r/27rT*) - ^'(1/2), (with r/2 instead of F for 
an s-wave and magnetic impurities^) , where ^'(a;) is the 
di-Gamma function and Tq is the pairing temperature 
in the absence of impurities. The ratio of the critical 
value of the scattering, : T*{Tcr) = 0, to Tq* is 
Tcr/T* = 7r/2eO-5^^2 ^ 9.88. 

The issue we address here both analytically and nu- 
merically is what is this ratio when u = 0(1) and ^ S> a, 
when the pairing problem involves incoherent fermions 
and near-gapless dynamical bosons and is very different 
from the d— wave version of the BCS theory. We find 
that in this situation the ratio Tcr/Tg is about 4 times 
larger than 0.88, i.e., the pairing is much less suppressed 
by impurities than in the weak coupling. This result is in 
agreement with the experiments which observed^i a sim- 
ilar reduction of the slope of T*{T) compared with the 
AG formula. 

Another issue that we consider here is how impurity 
scattering affects the angular dependence of the d— wave 
pairing gap. For a clean system, Abanov et al. have 
foundi that in the universal regime the form of the gap 
Ap{uj) is very close to cosp^; — cos py for all frequencies. 
We show that the cosp^ — cos py form holds in the pres- 
ence of impurities - the angular dependence only slightly 
changes with F. The implication is that both T* and the 
gap structure are robust towards impurities. 

The angular dependence of the gap, particularly in un- 
derdoped cuprates, has been the subject of intensive de- 
bates recently, and some ARPES data were interpreted 
as evidence for strong deviations from the cos Px — cos Py 
form. We emphasize in this regard that the position of 
the maximum of the spectral function Apitj) represents 
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the pairing gap Ap(w) only deep in the superconduct- 
ing state. At higher T, the position of the maximum in 
Apiuj) differs from Aj,(a;) because of damping induced 
by scattering off thermal bosons. In particular, even for 
a gap with a perfect cospx ~ cos Py form a maximum 
in Ap{uj = 0) is still present in some neighborhood of a 
node (a Fermi arc). In this regard, our result that the 
gap keeps cos Px — cos py form even in the presence of im- 
purities agrees with ARPES data by Campuzano et air- 
who detected this form at the lowest T. 

The effects of non-magnetic and magnetic impurities 
in superconductors with unconventional order parame- 
ters have been studied for high- Tc cuprateS )^°i^^'^^i^^'^'^i^^ 
non-cuprate superconductors, ^^^^'^'i^ and most recently 
for the pnictides.— For the cuprates, most of the stud- 
ies attributed a slow decrease of T* to the extended na- 
ture of the impurity potcntia l^^i^^'^"^ , but Monthoux and 
Pinesiii performed a numerical analysis of T* suppres- 
sion in YBCO by non-magnetic Ni impurities and found 
that the initial slope of T* is quite small even when im- 
purities are point-like scatterers. Very recently, Kemper 
et ali^ studied the effect of disorder using dynamical 
cluster approximation and quantum Monte Carlo, and 
found that ordinary pair-breaking by impurities is partly 
balanced by the impurity-induced enhancement of spin 
correlations which increases the pairing interaction me- 
diated by spin fluctuations 

Our result agree with Kemper et al>^ and also Graser 
et ali^ in that the origin of the flattening of T*(r) are 
magnetic strong-correlation effects. At the same time, 
we found that, in the universal regime, T* very weakly 
depends on the spin correlation length ^. In our theory, 
softness of T* suppression compared to weak-coupling 
AG theory is primarily associated with the strong fre- 
quency dependence of the pairing interaction. 

Theory. We follow earlier worki and consider 
fermions with a large FS and d— wave pairing mediated 
by ovcrdamped spin fluctuations. Wc add to earlier anal- 
ysis an isotropic, clastic scattering by point-like impu- 
rities. As customary for the pairing problem, wc in- 
troduce normal and anomalous Green's functions and 
self-energies and treat spin-fluctuation mediated pairing 



within the Eliashbcrg theory, by keeping self-energies but 
neglecting vertex corrections. For small and large m's, 
this approximation can be rigorously justified because 
vertex corrections are small in m or respectively. For 
u = 0(1), it can only be justified on the basis that vertex 
corrections are small numericallyi^ 

The set of equations includes fcrmionic and bosonic 
sclf-encrgies in the normal state, and the linearized equa- 
tion for the d— wave pairing vertex ^j^^{u}„i) (Ref[ll) 

Spj.(wm) == ttT* ^ J dp'fXpjZpr sign(w„w;„) , (1) 

) 



= -ttt* ^ / dp' r^^v^^ , , , 



,(3) 
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ia/0'+a^\Pf-p'f-Q\^ + \ALo\/n 



, (4) 



where = 3vf/{l6ua) and momenta pf in all formu- 
las are confined to the FS, because integration in the 
direction transverse to the FS has been carried out. In 
distinction to Ref. Qin (HD, © we also included impu- 
rity renormalization of , and of Matsubara energies 

i^m = TTT*{2m + 1), where F = {ui/irNf) sin^ S depends 
on impurity concentration rii, fermionic density of states 
Nf, and the impurity potential uq via tan 5 = iruoNf. 

The set of equations for <I> and S can be simplified in 
the usual way by introducing mass renormalization factor 
Zpf{iOm) and the pairing gap Apj.(a;m) via 



Zpj {iOm) 



Wm\ + r + Sp (w„i) 



\UJr, 



Pf\^mJ \ I \ _ ^Pfi'-^rii) 
, IXpj {UJrr,} ~ y (,, \ 



(5) 

Due to Aig symmetry of Ep^. and Big symmetry of 
the impurity renormalization of the pairing vertex van- 
ishes, i.e., $p^ = Using ([5]) we obtain from ^ 



ttT- 



Pf-Pf 



Z{uj^,Pf) 



AK„,p}) = 0. 



(6) 



We wrote the gap equation as an eigenvalue problem by 
moving all terms to one side and symmetrizing the kernel 
with respect to , ^^m' • 

Numerical solution. This linearized gap equation is 
solved numerically by presenting the FS integral as a 
sum, varying T and finding T* as the highest temper- 



ature where Eq. ([6]) is satisfied. The result is presented 
in Fig. [1] We clearly see a strong increase of the ratio 
Tct/Tq compared with a BCS d— wave superconductor. 
For u ^ 1, when Tq as a function of u has a maximum at 
about 0.02v//a, this ratio is nearly 4 times larger than 
in the BCS limit. This result is in a good quantitative 
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FIG. 1: (Color online) The onset temperature for the pair- 
ing, T* , vs. r for different values of u. Tq is the pairing 
temperature in the clean limit. Vf/a ~ leV in the cuprates, 
hence 10~^Vf/a ~ ImeV. We set £^/a = 4.8 for definiteness. 
The solid line is AG-type result for a BCS d— wave supercon- 
ductor, which in our case corresponds to the limit <^ 1. 
The key result in this figure is a progressive increase with u 
of the critical ratio Fcr/To , at which T* = 0. For u = 0(1), 
this ratio is about 4 times larger than in the BCS limit. The 
transition temperature is found with relative precision 10~^. 
The range of low T requires special care because for any fi- 
nite number of Matsubara points the curve T*{T) bends back 
towards the origin producing a spurious second solution. 




FIG. 2: (Color online) Left: the angular dependence of 
the quasiparticle renormalization factor Zpj{ijJm^o) and of 
the gap function Ap^(a;()) for it = 1, ^ = 4.8a (a) and 
u = 0.2, ^ = 9.8a (b). T*, F and tOm are in units of 10~''w//a. 
Observe that the angular dependences of the gap and Z (with 
subtracted constant shift due to F) changes very little be- 
tween F = 0(light lines) and F < Fcr(dark lines). Panel 
(c): frequency dependences of <^ = Ap^. (tJm) and Zpj{uirn)- 
Panel (d): the comparison of the numerically obtained fre- 
quency dependence of ^^^ois^m) with the analytic solution 
of Eq. (O in which momentum dependence of the gap near 
hot spot is neglected. 



agreement with the experiment in Ref. |8| and shows that 
the universal pairing scale in the cuprates is resistant to 
ordinary impurities. 

In Fig. [21 on the left we show the angular dependences 
of the quasiparticle renormalization factor Zp^^LOg) = 
1 + [Epj,(u;o) + r]/wo and of the gap function Apj:(wo) 
for clean and dirty cases, for luq = ttT* and different 
u and ^. In the clean case and u = 0(1) the angular 
dependence is quite close to cos 20 (or cosp^; — cospy). 
We see that the effect of the impurities on the angular 
dependence of the gap is quite small, i.e., cos 20 form 
is preserved in a dirty case. We verified that this holds 
for all F up to the critical value, and for all Matsubara 
frequencies. For completeness, in Fig. ^c) we show the 
frequency dependences of Zp, {ujm) and A 



{ujm), and in 



panel (d) we present $pj,(tJm) at 1^9 = which we will 
later compare with the analytical formula. 

Analytical reasoning. To understand the origin of the 
increase of Fcr/To* we analyze the equation for the pairing 
vertex ^pj{ujm) analytically at T ~ 0, i.e., we look for a 
solution near Tcr- To do this, we make an approximation 
and neglect angular dependence of the gap near a hot 
spot. The momentum integration along the FS then can 
be carried out analytically, and the equation for the gap 
at a hot spot becomes ID integral equation in frequency 
only. This equation is more easy to analyze than the orig- 
inal 3Z) integral equation. The approximation of the gap 
function by a constant near hot spots can be rigorously 
justified at small u (corrections are higher powers of u), 
but remains qualitatively valid up to u = 0(1) (Ref. [l|). 



The expression for Tq* in this approximation has been 
obtained earlier— - Tq w 0.13u(?j//a), with very weak 
dependence on ^ as long as > 1. For small w^, Tq is 
described by a BCS formula. 

Using ^pj+Q{ujjn) = -^pfi^^m), dropping the depen- 
dence of py near a hot spot, integrating over momentum 
in Eqs. (H]) - and rescahng variables, we obtain after 
some algebra the equation for ^{ujm) at T = in the 
form 



X 



dy- 



^4X^\y-x\ + l ^4A2|y + x-| + l 



,(7) 



where x,y are frequencies in units of a; = (3u/A)vf /a, 
r = Tcr/ijJ, and A = (2m^) (mass renormalization in the 
normal state is 1 -I- A). 

Weak coupling BCS limit corresponds to A ^ 1. In this 
limit, $(a;) becomes independent of x (and $ = A), the 
gap equation is solved in the same way as in AG theory, 
and the value of Vct/Tq « 0.88. We, however, are inter- 
ested in the opposite limit, when A > 1. We analyzed 
equation ([7]) by normalizing $(a;) to $(0) = 1, expand- 
ing at small x and at large x and extrapolating between 
the two limits. At small x, $ = 1 — 0(a;), at large x, 
$(a;) cx 1/^/x. We found that $(a;) is well approximated 
by ^{x) = 1/Vl + cx, where c is a constant which de- 
pends on A and F. The error is less than 2% for all x 
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and for all A which we considered. We also checked that 
that this solution is not an artefact: if we use this ^{x) 
as an input and run iterations, ^{x) rapidly converges. 
The critical value of F is then obtained by substituting 
this form back into ([7]) and solving for $(0) = 1. 

Carrying out this procedure, we found that c increases 
with increasing A = u^, i.e., when the correlation length 
increases and the pairing problem becomes more and 
more non-BCS, the gap function gets confined to pro- 
gressively smaller frequencies. The c(A) increases from 
c(l) = 0.66, through c(2) = 1.04 and c(5) = 1.7, to 
c{oo) = 2.31. In Fig. [2Jd) we compare our approximate 
analytic ^{uJm) for u = 0.2, ^ = 9.8a (A w 2) with the 
numerical (f>pj. (wm). We see that the agreement is ex- 
pectedly not prefect, but generic trends of the frequency 
dependence is captured by the approximate solution. 

Substituting $(x) = 1/^/1 + cx back into ([7]) and solv- 
ing for $(0) = 1, we found that F progressively increases 
as A gets larger, from X~^e~^^'^ at small A to 0.3 for A = 1 
and to 0.46 for A = oo. 

The monotonic increase of the value of F with increas- 
ing A is a tricky effect. One could expect that the confine- 
ment of to smaller x as A increases and the increase 
of the self-energy tend to reduce F simply because typi- 
cal frequencies get smaller. However, as A increases, the 
interaction strength also increases, and this tends to in- 
crease F because F appears in the denominator in the in- 
tegral for $(0) = 1, and larger F arc required to balance 
the increase of the interaction. We compared the two 
effects and found that increase of the interaction over- 
shadows other effects and is the origin of the growth of 
F with increasing A. 

We next compared the growth of Tcr ~ a;r(A) and 
the grown of Tq . The latter also scales as ui with A- 
depcndcnt prefactor (Ref. [23[ ). This prefactor increases 
with increasing A, but its A-depcndencc is very weak: it 
changes by less than 5% between A = 1 and A = oo. As 



a result, the A-dependence of the ratio Tct/Tq predom- 
inantly comes from F(A), which, we remind, increases 
with A. Inserting the numbers, we find that the ratio 
Tcr/T* becomes 2.0 for A = 1; 2.37 for A = 2; 2.47 for 
A = 5, and 2.74 for A = oo. The scale of the increase is 
quite consistent with what we found numerically in Fig. [T] 
by solving the full 3D integral equation in momentum and 
frequency. 

Summary. In this paper we considered the effect 
of non-magnetic impurities on the onset temperature T* 
for the d— wave pairing in spin-fluctuation scenario for 
the cuprates. Non-magnetic impurities are pair-breaking 
for c?— wave superconductivity, and one should expect a 
reduction of T* due to impurities. In weak-coupling, T* 
falls off rapidly, following the AG curve. 

We analyzed the effect of impurities in the intermedi- 
ate coupling regime when the magnetic correlation length 
^/a > 1, the dimensionless coupling u is 0(1), and the 
pairing problem is qualitatively different from BCS. In 
the clean limit, T* in this parameter range weakly de- 
pends on ^ and u and is approximately 0.02vf/a. We 
found that this universal pairing scale is quite robust with 
respect to impurities: the critical value of the scattering 
rate Tcr needed to bring T* down to zero is about 4 times 
larger than in the weak coupling. This implies that the 
slope of the initial reduction of T* is weaker by about the 
same factor than in the weak coupling. This reduction 
of the slope agrees with the experiments^ and with ear- 
lier work by Monthoux and Pincsiii on T* suppression 
in YBCO due to non-magnetic Ni impurities. Wc also 
analyzed the angular dependence of the gap and found 
that it is little affected by impurities. 
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Three of us (A.B.V., M.G.V. and A.V.Ch.) arc thankful 
to the Aspen Center for Physics for hospitality. 
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